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ABSTRACT 



As part of a program to understand disk dispersal and the interplay be- 
tween circumstellar disks and X-ray emission, we present new high-resolution 
mid-infrared imaging, high-resolution optical spectroscopy, and Chandra grating 
X-ray spectroscopy of the weak-lined T Tauri star DoAr 21. DoAr 21 (age < 10 6 
yr and mass ~ 2.2 M Q based on evolutionary tracks) is a strong X-ray emitter, 
with conflicting evidence in the literature about its disk properties. It shows 
weak but broad Ha emission (reported here for the first time since the 1950s); 
polarimetric variability; PAH and H 2 emission; and a strong, spatially-resolved 
24-/im excess in archival Spitzer photometry. Gemini sub-arcsecond-resolution 
9-18 /im images show that there is little or no excess mid-infrared emission within 
100 AU of the star; the excess emission is extended over several arcseconds and 
is quite asymmetric. The extended emission is bright in the UV-excited A = 11.3 
/im PAH emission feature. A new high-resolution X-ray grating spectrum from 
Chandra shows that the stellar X-ray emission is very hard and dominated by 
continuum emission; it is well-fit by a multi-temperature thermal model, typical 
of hard coronal sources, and shows no evidence of unusually high densities. A 
flare during the X-ray observation shows a temperature approaching 10 s K. We 
argue that the far-ultraviolet emission from the transition region is sufficient to 
excite the observed extended PAH and continuum emission, and that the H 2 
emission may be similarly extended and excited. While this extended emission 
may be a disk in the final stages of clearing, it also could be more akin to a small- 
scale photodissociation region than a protoplanetary disk, highlighting both the 
very young ages (< 10 6 yr) at which some stars are found without disks, and the 
extreme radiation environment around even late-type pre-main-sequence stars. 
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Introduction 



An open question in early stellar evolution and planet formation is the process and 
timescale by which stars disperse their protoplanetary disks. Central to this question is the 
lifetime of the gas component of circumstellar disks. The gas dominates the disk mass, but 
since small dust grains dominate the infrared opacity, dispersal or growth of these grains can 
significantly reduce or eliminate the excess infrared emission that is often taken to be the 
primary signature of a circumstellar disk, even if substantial gas is still present. 

Related to the question of gas disk survival and evolution is the influence of X-ray 
emission on the circumstellar gas, and vice versa. T Tauri stars , low-mass pre-main-sequence 



stars, are often strong X-ray emitters (see, e.g., the review of lFeigelson fc Montmerld 11999 



and references therein). A detailed understanding of the X-ray emission from these young 
stars is important both for an understanding of the early evolution of the stars themselves, 
and due to the X-rays' impact on the surrounding environment. T Tauri stars are often 
surrounded by disks of gas and dust, from which planets presumably form later in the pre- 
main-sequence phase. The physical state of these disks is strongly influence d by the star's X- 
ray e mission. The X-rays are the dominant source of ionization of the disk (jlgea fc Glassgold 



19991 ). they strongly influence the di sk chemistry (jMalonev et al 



1996h . and they can be 



important in disk photoevaporation (jGorti fc Hollenbach 2009 : Drake et al. 2009 ). X-ray 



flares from the young Sun may have p roduced the isotopic anomalies seen today in meteorites 
( Feigelson. Garmire. &: Pravdo!l20o"2l ). 



Similarly, the presence of circumstellar gas may play a role in shaping or modifying the 
star's X-ray emission. The model for the X-ray emission of the vast majority of T Tauri stars 
is that the X-rays are produced by solar-type magnetic activity scaled up by rapid rotation 



and / or longer convective turnover times (e.g., iFeigelson fc Montmerld Il999l ; iPreibisch et al. 



20051 ). In a small number of systems, there is also evidence that the X-rays are influenced by 
i nteraction between a stellar magnetic field and a circumstellar di sk, perhaps via accretion 
( jKastner et al.l 120021 : IStelzer fc Schmittl I2004J ; ISchmitt et al.l 120051 ) . The unprecedented X- 
ray spectral resolution and sensitivity of the Chandra X-ray satellite offers an opportunity 
to place much better constraints on the temperature and density of the X-ray emitting gas, 
thereby providing potentially important information about the X-ray emission mechanism(s), 
and about the X-ray radiation field incident on the circumstellar disk. 
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To date, however, only a modest number of high-spectral-resolution X-ray observations 
of pre-main-sequence stars have been published; even the distances to the nearest star- 
forming regions result in X-ray fluxes that make high-resolution o bservations prohibit ively 
long. Ther e are a handful of not able exceptions, however. TW Hya (IKastner et al.ll2002l ) and 
HD 98800 (IKastner et al.ll2004j ) both lie in t he nearby TW Hya ass ociation, with a distance 
of roughly 50 pc and an age of 5-15 Myr (IWeintraub et al.ll2000l . and references therein). 
The Chandra X-ray spectrum of TW Hya yielded a surprisingly low temperature for the 
X-ray e mitting gas, and sug gestions of very high densities in the X-ray emitting regions, 
leading IKastner et al.l (120021 ) to suggest that the X-ray emission is related to accretion from 
the star's circumstellar disk. IStefzer fc Schmittl (12004 ) reached a similar conclusion based 
on density diagnostics in TW Hya's XMM-Newton spectrum, and also attributed the low Fe 
abundance seen in the X-ra y emitting gas t o dep letion of some elements from the gas phase 
onto dust grains in the disk. ISchmitt et al.l (120051 ) presented similar evidence for high-density 
X-ray emitting gas around the classical T Tauri star BP Tau. 



In contrast, the X-ra y spectra of both HP 98800 (IKastner et al.ll2004j ) and the similarly- 
aged (~ 10 Myr) PZ Tel (lArgirofh et al.ll2004l ) are quite similar to spectra of older solar-type 
stars and show no evidence for unusually high densities, suggesting that their X-ray activity 
i s, to first order, scaled-up solar activi ty. Solar-like X-ray activity is also seen in AB Dor 
(jSanz-Forcada. Maggio. fe Micelal 120031 ). which is near or recently arrived on the zero-age 
main sequence. Notably, TW Hya and BP Tau show clear evidence of active accretion from 
their circumstellar disks, while the stars with more solar-like X-ray emission do not have 
disksE In general, while accretion may modify the X-ray spectrum in some young stars, 
studies of the X-ray emission fro m large samples of young stars show that accretion is not a 
primary driver of X-ray activity (jPreibisch et al.ll2005l ; IStassun et al.ll2006l ). 

In the context of understanding both disk evolution and the interplay between disks and 
X-rays, the weak-lined T Tauri star DoAr 21 (a.k.a. V2246 Oph, Elias 2-14, VSSG 23, GSS 
23, ROXs 8, Haro 1-6, YLW 26, HBC 637, ROXR1 13) presents an interesting case study. It 
has no strong evidence for accretion or a massive disk, and yet it is one of the few T Tauri 
stars around which emission f rom H ? and polycy clic aromatic hydrocarbon (PAH) molecules 



has been detected (Sec. I2.2p . iBary et al.l (120021 ) suggest that it may r etain a gas-only disk 
in which all of the dust has coagulated into larger grains; Cieza et al. ( 20081 ) classify it as a 
"transition disk," one in the process of being cleared from the inside out. In addition, it is 
the most luminous X-ray source in p Oph, making it an excellent laboratory to study the 



1 Only the southern comp onent ( HD 98800 B) of the wide binary pair in HD 98800 has circumstellar 
material ( Koerner et al. 2000f ). while Kastner et al. ( 2004 ) show clearly that the X-ray spectrum they detect 
and analyze is from the diskless northern component, HD 98800 A. 
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influence of X-rays and circumstellar material. 

Here we present new Gemini mid-infrared images, high-resolution optical spectra, and 
Chandra high- resolution grating spectra of DoAr 21. The primary contrast of DoAr 21 with 
the handful of pre-main-sequence stars previously observed at high X-ray spectral resolution 
(Fig. [1]) is the combination of its youth and lack of a strong accretion signature. DoAr 21 is 



embedded in the p Oph star- forming region (see IWilking. Gagne. fc Allen! 120081 for a recent 
review), whereas a number of the others (e.g., TW Hya, Hen 3-600, MP Mus, and V4046 
Sgr, with ages of ~ 10-30 Myr) no longer lie in regions of active star formation, with all 
surrounding molecular gas having been dissipated. With an estimated age of < 1 Myr, DoAr 
21 is younger than these other stars. A few of the stars observed with high X-ray spectral 
resolution (e.g. T Tau, BP Tau, RU Lup) are of a similar age to DoAr 21, but all are classical 
T Tauri stars with strong accretion signatures. The central questions here, then, are how 
much difference the young age and lack of strong accretion make in the X-ray emission of 
DoAr 21; what the distribution and physical conditions of its circumstellar material are; and 
how the high-energy radiation field incident on the circumstellar material may influence its 
properties. 

Below we discuss the stellar and circumstellar properties of DoAr 21 (Sec. [2]), showing 
new images of its circumstellar material, which, surprisingly, is not present near the star 
but which extends to several hundred AU. We present new high-resolution optical spectra 
which show that DoAr 21 is probably not accreting, but which are consistent with strong 
chromospheric activity. In this light, we then present our new X-ray observations (Sec. [3]) 
and analyze the X-ray emission in detail through emission-line-ratio fits, modeling of the 
entire spectrum, and analysis of time variability of the emission before, during, and after a 
large flare. In Sec. HI we discuss the influence of the star's X-ray and ultraviolet emission 
on its circumstellar environment. Based on the level of activity suggested by our X-ray 
observations, we estimate the far-ultraviolet (FUV) flux, suggesting that it is sufficient to 
excite the observed PAH emission even at hundreds of AU from the star. Finally, we present 
a possible geometry for the H2 emission, suggesting that it may be more akin to a photo- 
dissociation region (PDR) than a circumstellar disk. 



2. DoAr 21's stellar properties and circumstellar environment 

Before considering our new X-ray data for DoAr 21, here we consider new and existing 
data from other wavelengths, including new mid-infrared imaging and high-resolution optical 
spectroscopy. We also consider infrared and optical photometry and spectroscopy from 
the literature in order to build up a more complete picture of the stellar properties and 
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Fig. 1. — HR diagram showing the position of DoAr 21 (star symbols) compared to other 
pre-main- sequence stars that have been observed with the Chandra or XMM-Newton grat- 
ings. The dotted star shows the luminosity of each star in DoAr 21 if it is an equal-mass 
binary. Filled symbols are classical T Tauri stars, and open symbols are weak-lined T 
Tauri stars. DoAr 21 is both younger (~ 4 x 10 5 yr) and more massive (~ 2.2 M ) 
than most of the T T auri stars observed preyiously . The pre- main-sequence e v olutio nary 



tracks shown are fromlD'Antona fc Mazzitellil (119971 ) . Both the iPalla fc Stahlerl (119991 ) and 



Siess. Dufour. Sz Forestini (boooh tracks result in slightly older absolute ages in general, but 



similar relative ages. 



- 6- 



circumstellar environment of DoAr 21. 



2.1. Stellar properties and multiplicity 



Spectral types in the literatur e are KO (Bouvier &: Appenzellerlll992l ; Martin et al.lll99£ 
from o ptical spectra and K0-K 2 (jLuhman fc Riekd Il999f ) from near-infrared spectra. Fol- 
lowing 



of 5080 K 



Luhman &: Riekd (Il999l ). we adopt a spectral type of Kl and effective temperature 



T he dis tance to DoAr 21 from a VLBA parallax measurement is 122 ± 6 pc 



(ILoinard et al.l 120081 1 . Using this distance and our fit to the de-reddened spectral energy 
distribution (Sec. 12.21) . we estimate the stellar luminosity to be 11.7 ± 0.8 L , where the 
uncertainty arises primarily from the large photometric variability and somewhat uncertain 
extinction. 



The VLBA observations of ILoinard et al.l (120081 ) show DoAr 21 to be a binary, with a 
projected separation of 5 mas (0.6 AU) in one of their observations. The orbital parameters 
are uncertain, but Loinard et al. estimate the semimajor axis to be 1-2 AU. Since the binary 
pair has been detected only at radio wavelengths thus far, there is no information about the 
mass ratio or the optical luminosity ratio. Ascribing all of the luminosity to the primary 
star (solid-lined star in Fig. [TJ yields a mass of ~ 2.2 M & and an age of ~ 4 x 10 5 yr. If the 
luminosity is split equally between the two stars (dotted star in Fig. [TJ), the mass of each star 
is ~ 1.8 M & and the age is ~ 8 x 10 5 yr. There is some evidence in the optical spectra (Sec. 
12. 4p that the secondary is contributing to the observed line width, suggesting that its optical 
luminosity is non- negligible; thus the latter set of parameters may be closer to reality. We 
note that DoAr 21 is relatively massive among T Tauri stars, and it has a K spectral type 
only because of its very young age. By an age of 3-4 Myr, it will move into the spectral type 
range of Herbig Ae/Be stars, and it will arrive on the main sequence as a late B or early A 
star. 



The spectral type of B2 found in some references, and as of this writing in Simbad, is from lChini et all 



(|1977I ). and is based on B, R, I, and K photometric colors. The authors did not find a unique combination 
of spectral type and extinction that fit the data; assuming an extinction of Ay = 10.1 they found a s pectral 
type o f B2, and assuming Ay = 6.6-6.2 they found a spectral type of G5-K0. A follow-up paper by IChini 
(|198ll ) used a broader range of colors and found Ay — 6.2 and a spectral type of K0, consistent with what 
we find here. 
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2.2. Spectral energy distribution and infrared excess 



The evidence in the literature for circumstellar material around DoAr 21 (and accretion 



of su c h material) is mixed, and DoAr 2 1 has often been consid ered (e.g., lLada &: Wilking 
19841 : lAndre et all Il992t IPreibischl Il999l : iBontemps et all bpOlh to be a Class III or disk- 
less T Tauri star, though s ome sources have stated that it has an infrared excess (e.g., 
Bouvier fc Appenzellerlll992l ). and it has recently been classified as a "transition di sk" sys- 



tem, with a cleared inner disk but optically thick outer disk (e.g., ICieza et al.ll2008l ). 



In light of the lack of a strong infrared excess, t he detection of em i ssion features from 
PAHs ( banner. Brooke, fc Tokunagal Il995l ) and H 2 fearv et all 1200.4 bitner et all booah 
was surprising, raising the possibility that a substantial gas component of th e disk is sti ll 
present. Early observations showe d Ha with bright and var i able em i ssion (jHarol Il949l ). 
and then later weakly in emission (IDolidze fc Arakelyanl Il959l ; iHidaiatl Il96lh . Since then . 



however, no published observations of the star have shown Ha emission ( Rvdgren et al 



Montmerle et al. 



1992: 



19831: 



Martin et al 



Feigelson & Montmerldll985l : IWilking et al.lll987 



1976 



Bouvier &: Appenze 



1998h . Br7 is in absorption with no near- infrared veiling ( Luhman fc Rieke 



ler 



19991 ). Our new high-resolution optical spectra (Sec. I2.4[) show that the photospheric Ha 
absorption line is partly filled in with emission and is surrounded by weak but very broad 
(~ 300 km s _1 FWHM) emission wings, revealing a weak and variable emission component 
that is consistent with the Ha emission seen in older, coronally-active low-mass stars and 
which thus suggests that there is little to no accretion onto DoAr 21. 

The spectral energy distribution (SED) of DoAr 21 gives some clues to the nature 
of its circumstellar environment. We searched the literature for available photometry of 
DoAr 21, and assembled the SED from A = 0.39 /im (U band) to A = 1.3 mm. The 
observed broad-band colors are much redder than photospheric colors for a Kl star. To 
de-redden the photometry and to determine whether or not DoAr 21 shows any evidence 
of infrared excess, a signature of the presence of a circ umstellar or circumbin ary disk, we 
compared the SED of DoAr 21 to a model photos phere (IBuser fc Kuruc z 1992) with T c g = 
5000 K (Fig. HP - We adopt t he extinction law of Weingartner fc Draine (j2001a ). with the 
modifications of Draine J2003h . This extinction law has been shown to be a better fit to mid- 
infra r ed data than a power - law extrapolati on from the near-infrared (jRoman-Zuniga et al. 

ir 



20071 : iChapman et all 120081 : iMcClure 

Rv = Ay/E{B — V) = 4.2, determined from the relationship Ry 



20091 ) . We use a ratio of total-to-selective extinction 

5.7A max (IVrba et al. 



19931). where Ar 



0.74 (Urn is the wavelength of maximum observed polarization for DoAr 



21 (IMartin et al.lll992l). This value o f Ry is consistent with results found for other lines of 
sight in Ophiuchus (IVrba et al.lll993l ) and yields a notably better fit to the UBVRI data 
than the standard interstellar value of Ry = 3.1. Physically, a larger Ry corresponds to 
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larger grain sizes, as expected in a dense cloud. Using this extinction law, we de-reddened 
the photometry to find the best fit of the model photosphere to the V, R, and / data, 
as the U and B data show significant variability and the longer-wavelength data may be 
affected by any infrared excess. This procedure yielded Ay = 6.2 mag, and combined 
with our adopted distance, a photospheric luminosity of 11.7 ± 0.8 L & . This corresponds 
to a stellar radius of 4.6 ± 0.2 R & if all of the luminosity is attributed to one star. The 
adopted Ay corresponds to Aj = 1.8. We note that this is consistent with the neutral 
hydrogen column density of pa 10 22 atoms cm -2 found from the X-ray spectrum (Sec. 
3]); from a s a mple of 20 sources with good X-ray and infrared data in the p Oph cloud, 
Vuong et al.l (120031 ) find N H /Aj = 5.6 ± 0.4 x 10 21 cm -2 mag -1 . Given the differences 



in extinction law s used, our adopt e d ext inction is reasonably con sistent with that found 
by oth er authors; iLuhman fc Riekd (119991 ) find Aj = 1.6 mag, and iBouvier &: Appenzeller 
(119921 ) find Ay = 6.6. The de-redde ned photometric d ata ar e plotted in Figur e [2J the 
millimeter-wavelength upper limits of I Andre et al.l (jl990j) and ICieza et al.l (120081 ) are not 
shown. 



While the observed optical and near-infrared variability (see Section 14.1.51) makes it 
impossible to fit every data point, the photometry is broadly consistent with the model 
photosphere within the uncertainties for A < 5 /im. There is some evidence for excess 
emission at U and B, which may be due to flaring as discussed below. While the brightest U 
and B points are clearly inconsistent with photospheric emission, the fainter points lie close 
to the model photosphere, consistent with a scenario in which DoAr 21 is flaring much of 
the time. 

At longer wavelengths, there is a clear infrared excess at A = 7-24 //m. These points 
are quite insensitive to extinction corrections, since the extinction is so small, and to the 
determination of the star's effective temperature, since they lie on the Rayleigh- Jeans tail 
of the photospheric emission. The space-based infrared data are highlighted in the figure as 
squares (orange in the on-line edition), as they are not affected by the Earth's atmosphere 
and are m easured on a uni f orm p hotometric system; the A = 6.7 and 14 /zm points are ISO 
data from iBontemps et al.l (12001 ), whil e the other A = 3.6-24 fim squares are Spitzer IRAC 
and MIPS data from lCieza et al.l (120071 ). While the presence of an excess shortward of 7 jum 
is debatable, the 8-24 /im points are well in excess of the photospheric emission, indicating 
the presence of some circumstellar material within a radius of 4"-9", the aperture size of 
the mid-infrared observations. This circumstellar material apparently does not have a large 
mass in small du st grains, however, given the very sensitive 3a flux upper limits of < 4 mJy 
at A = 1300 fan (Undre et al.lll990h and < 18 mJy at A = 850 /an Jcieza et al.ll2008h . 
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Fig. 2. — The spectral energy distribution of DoAr 21, from ultraviolet through infrared 
wavelengths. Spitzer and ISO data are shown as filled squares (orange in the on-line edition), 
and the 850 /im and 1.3 mm upper limits are not shown. The photometry has been de- 
reddened with Ay = 6.2 and Ry = 4.2 as described in the text. Overplotted is a model 
photosphere (solid line); the star shows a clear infrared excess at A > 8 /xm. 
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2.3. Mid-infrared imaging 

2.3.1. Archival Spitzer images 
DoAr 21 was observed by the Spitzer Space Telescope with both IRAC (A = 3.6, 4.5, 5.8, 



8.0 jum) and MIPS (A = 2 4, 70, and 160 /im) as part of the c2d Legacy program (jCieza et al. 



20071; iPadgett et all 20081 ). PSF-fitting photometry of DoAr 21 from 3.6-24 /im is reported 
by lCieza et al.1 (120071 ) and plotted in Figure [2j here we concentrate on the images themselves. 
Figure [3] shows the A = 24 /im image of the region around DoAr 21; DoAr 21 itself is clearly 
detected. Visual examination of the A = 70 /im images shows a source at the position of 
DoAr 21, but the noise is too high for accurate photometry. 

Surprisingly, the 24-/xm emission from DoAr 21 is somewhat extended. Close inspection 
of the 24-/im MIPS images reveals that the profile of DoAr 21 is different from that of other 
stars in the field. Since Spitzer is diffraction-limited at 24 /an, the point-spread-function 
(PSF) is a modified Airy pattern, and indeed many of the sources in Figure [3] (including 
sources that are fainter than DoAr 21 and/or projected on comparable or brighter background 
emission) show the presence of the first Airy ring. DoAr 21 does not, however, indicating 
that it is not a point source at 24 /im. To check this, we compared the azimuthally-averaged 
radial profile of DoAr 21 to that of a star in the same image (GSS 26), and also to a 
sample empirical on-orbit PSF (constructed from many stars in a MIPS mosaic image of the 
Trapezium) supplied by the Spitzer Science Center (SSC). DoAr 21 is clearly more extended 
than the other stars, both of which show the signature of the first Airy ring at radii of 3-5 
pixels, while DoAr 21 does not. (The diffraction limit of Spitzer at 24 /im is 7'.'4, while the 
MIPS mosaic pixel size is 2'.'45, so the first null occurs at a radius of 3 pixels.) We examined 
three different MIPS mosaic images, and DoAr 21 appears extended in all of them. We also 
examined the individual baseline-calibrated data (BCD) frames from the SSC archive (from 
which the mosaics were constructed), and found that DoAr 21 appears extended in each 
frame, regardless of its position on the chip, while other stars in the frame do not. Thus, 
the broader profile of DoAr 21 is not an artifact of the mosaicing process. 

To estimate the radial extent of the emission, we convolved an oversampled MIPS PSF 
with Gaussians of different radii, resampled them at the MIPS pixel scale, and compared 
their radial profiles to that of DoAr 21. The profile of DoAr 21 is roughly consistent with 
a Gaussian source with a sigma of 0.75-1 MIPS pixel (l'!8-2!!5). At the 122 pc distance of 
DoAr 21, this corresponds to a FWHM of 500-700 AU. Thus, the 24-/tm emission appears 
to be extended to radii of 250-350 AU from the star. 
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Fig. 3. — Spitzer MIPS 24-/zm image of the field around DoAr 21; DoAr 21 is clearly 
detected, indicating a substantial amount of circumstellar material. Unlike other stellar 
sources in the field, DoAr 21 does not show the first Airy ring, indicating that the 24-/zm 
emission is marginally resolved. The B2 star HD 147889 (lower right) excites the emission 
from the bright ridge to the northeast of DoAr 21. 
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2.3.2. New Gemini images 

Motivated by the extended emission seen in the Spitzer images, we observed DoAr 21 
with the T-ReCS camera on Gemini South in June and July 2007. The observations were 
made through filters with central wavelengths (and 50% transmission widths) in /im of 8.6 
(AA = 0.43), 10.4 (AA = 1.02), 11.3 (AA = 0.61), and 18.3 (AA = 1.51). The 8.6 and 
11.3 /im filters are centered on PAH features, while the 10.4 /im Si-4 and 18.3 /im Qa filters 
largely measure continuum emission. In particular, the Si-4 filter acts as a nearby continuum 
reference for any PAH emission. 

The data were taken with the standard chop and nod technique to remove sky back- 
ground emission, and were processed with the standard Gemini IRAF pipelined The chop 
throw was 15" at a position angle of 130°, parallel to the ridge of 24-/im emission seen in 
the Spitzer images (Fig. [3]); this chop direction was chosen to minimize variation of the 
background emission. 

The data were flux-calibrated with observations of HD 136422 (K5 III) at a similar 
airmass and near in time to the DoAr 21 observations in each filter. Fluxes for HD 136422 



were taken from the calibrated template spectrum created by ICohen et al.l (119991 ) . Because 
the target and calibrator were observed at similar, low (< 1.3) airmasses and near in time 
to each other, we did not attempt to correct for differential atmospheric extinction between 
target and calibrator observations. 

The spatial resolution of our observations, as measured by the FWHM of standard stars 
observed at similar airmass to DoAr 21, was in the 8.6-11.3 /im observations, and 

0'/58 at 18.3 fim. This is near the diffraction limit of Gemini South, and is well-sampled by 
the T-ReCS pixel scale of 0"09 pixeD 1 . 

The resulting images are shown in Figure HI Irregular, extended emission is seen at 
all wavelengths. The emission is brighter in the filters containing PAH features (8.6 and 
especially 11.3 fim) than in the 10.4 /im continuum, suggesting that the extended material 
shows PAH emission. The brightest off-source emission is a spot roughly l'/l north of DoAr 
21, and there is a ridge or partial arc of emission to the north and west, seen in all filters. 

The central stellar source in the 8.6-11.3 /im DoAr 21 images appears point-like at our 
resolution; it has a FWHM similar to that of the photometric calibrator in the corresponding 
band. The slight extension of the PSF at PA ~ 230° in the 8.6 /im images is also seen in 



3 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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Fig. 4. — Gemini T-ReCS images of DoAr 21 at four different wavelengths. The 8.6- and 
11.3-/zm bands are centered on PAH emission features, while the other two bands trace 
continuum emission. The grayscale shows the emission displayed with a logarithmic stretch; 
the 18.3-/zm image is smoothed with a 3-pixel FWHM Gaussian. The circle in the lower left 
of each image shows the FWHM of the PSF measured from a comparison star observed at 
similar airmass. 
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the calibrator. The central source in the 18.3 /im image has FWHM of O'.'6-l'.'O; given the 
much lower signal-to-noise ratio at this wavelength, it is difficult to tell if the central source 
is extended, though certainly there is significant extended emission beyond this. 

The extended emission is quite asymmetric, unlike the emission expected from a disk. 
There is a partial arc or ring to the north and west (seen most prominently at 11.3 and 18.3 
/im) at roughly l'/l from the star. There are bright knots at position angles of roughly 0° 
and 140°, with the northern knot appearing both brighter and more concentrated, and the 
SE knot being visible to larger radii, especially at 18.3 /im. As shown in the star-subtracted 
azimuthal profile in Figure [5], there is a roughly 4:1 contrast between the brightest and 
faintest position angles. 

Archival HST WFPC/2 images of DoAr 21 in the F606W and F814W filters show a 
faint filament extending 8"-10" southeast from the star; the filament is not apparent in 
im ages in the F104 2 M fil ter. There is no obvious brightness gradient along the filament, 



so 



Stapelfeldt et al.l (120091 ) suggest that it is foreground or background material illuminated 



by the star, rather than a jet associ ated with the star. This filament is also visible in the 



narrow-band near- infrared images of iGomez et al.l (120031 ). in filters both on and off the H 2 



1-0 S(l) line. Comparison of the WFPC/2, near-IR, and Gemini images shows that this 
filament overlaps the extension seen in Figure H] at PA 140°. Thus, the circumstellar material 
apparently extends to even larger radii than seen in our images, at least in some directions. 
We discuss this point further in Section HI 

To quantify the spatial distribution of circumstellar material, and in particular to see 
whether there is any infrared excess emission that is not spatially resolved, we measured the 
flux in our images in apertures of different radii. We chose aperture radii of 0'/82, 2? 75, and 
6". The smallest radius corresponds to a physical radius of 100 AU, a typical disk radius 



for a young star; the intermediate radius matches the aperture used by lHanner et al.l (119951 ) 
for their mid-infrared spectroscopy; and the largest radius encompasses all the detected flux 
above the background level in our images, probing roughly the same spatial scale as the 
MIPS 24-/im images. These fluxes are quoted in Table [1] and plotted in Figure [61 There is 
little infrared excess within 100 AU of the star, while the photometry of the more extended 
emission is roughly consistent with previous photometry. 



Given the clear presence of a significant amount of circumstellar material, but evidence 
that little of it is close to the star, we now turn to the question of whether there is any 
evidence for accretion onto DoAr 21. 
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Fig. 5. — Azimuthal profile of the 11.3-//m emission between radii of 100 AU and 550 AU, 
showing the 4:1 contrast in brightness at different position angles. The profiles in the 8.6 
/im and 10.4 /im images show similar structure and contrast, though at lower flux levels. 
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Fig. 6. — A close-up of the mid- infrared SED, with the model photosphere and de-reddened 
photometry as in Figure El The broad-band Spitzer an d ISO data ar e show n as filled squares 
(orange in the on-line edition), and the spectrum from lHanner et al.l (I1995I ) is shown as open 
squares (green in the on-line edition), showing an emission feature at 11.3 /mi attributed 
to PAHs. De-reddened photometric measurements from new Gemini narrow-band images 
(wavelengths 8.6, 10.4, and 11.3 /im) are overplotted, with horizontal bars showing the 
bandwidth of the filters used and vertical error bars showing the assumed 10% uncertainty. 
The open triangles (blue in the on-line edition) show measurements within 100 AU (0'.'82) 
of the star, showing little or no excess in the continuum at 10.4 /im. The open diamonds 
(purple in th e on-line edit i on) sh ow measurements in a 2?75 radius aperture, the same as 
that used by lHanner et al.l (119951 ). 
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2.4. Optical spectrum 



As noted in Section 12.21 the Ha emission from DoAr 21 has shown varying behavior 
at different times in the past, and thus additional spectroscopic observations are useful to 
probe this behavior. In addition, DoAr 21 's stellar rotation p eriod is unknown, as it did 
not show periodic photometric variability in the observations of iBouvier et all (119881 ); thus, 
measurement of its projected rotational velocity v sini is useful. 

To probe the Ha emission and stellar rotation of DoAr 21, we observed it with the echelle 
spectrograph on the CTIO Blanco 4-meter telescope on 2002 June 12 and 14, and on the 
four consecutive nights 2003 April 11-14. The spectrograph setup gave wavelength coverage 
of approximately 5,000-8,000 A with 0.08 A/pixel. With a 1" slit width, this gave a spectral 
resolving power (as measured from the FWHM of narrow lines in a ThAr comparison lamp 
spectrum) of R = 40, 000 at the Ha line. Exposure times were 20 minutes on 2002 June 
12, 15 minutes on 2003 April 11 and 30 minutes on the other four nights. The data were 
reduced using standard routines for echelle spectra in IRAF. 

Our spectra show broad absorption lines. There is no obvious evidence in the spectra 
of doubled lines from a binary companion, though the broad lines and relatively low signal- 
to-noise ratio limit our ability to detect other lines. The Li I 6708 A a bsorption line has a n 
equivalent width of 320 ± 20 mA, consistent with the measurement of iMartm et all ( 19981 ). 
From our spectra we measure a rad ial velocity of —6 ±4 km s _1 , consistent with the velocity 
of —4.6 ± 3.3 km s _1 measu red by Massarotti et all (120051 ). and with the p Oph young-star 
mean of —6.3 ± 1.5 km s _1 (|Pratoll2007l ). Comparison of our spectra with a Kl V standard 
spectrum artificially broadened to various velocities gives a projected rotational velocity 
vsini = 80 ± 10 km s _1 for DoAr 21. Combined with the stellar radius of 4.6 R Q estimated 
from its effective temperature and luminosity (attributing all the luminosity to one star), 



Table 1. Photometry of DoAr 21 


Filter 


Flux within 0"82 


Flux within 2"75 


Flux within 6" 




(mJy) 


(mJy) 


(mJy) 


8.6 /iin (PAH) 


313 


525 


596 


10.4 /iin (Si-4) 


176 


300 


347 


11.3 //m (PAH) 


234 


680 


816 


18.3 /im (Qa) 


165 


607 


728 
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this gives a rotational period of around 3 days, or less if the star is not viewed equator-onQ 



Given the discovery that DoAr 21 is a bina ry (Sec. 12. ip, how e ver, t here is an alternate 
interpretation of the width of its spectral lines. iMassarotti et ajj (120051 ) report vsini = 29 
km s _1 for DoAr 21. The origin of the discrepancy between their determination and ours 
is unclear, but it could result from observing a binary at different orbital phases, so that 
sometimes the lines are broa dened (though no t clear ly doubled) due to different radial ve- 
locities of the two stars. The IMassarotti et al.l (120051 ) spectra, measured at similar spectral 
resolution but with less spectral coverage and lower S/N, show some variation in the width 
of the cross-correlation peak, which is a measure of the line width, occasionally showing 
broader peaks (D. Latham, personal communication). Our spectra show consistently broad 
lines, and none of our six spectra is con s istent with v sini as low as 29 km s _1 . We repeated 
the procedure used by IMassarotti et al.l (120051 ) on our spectra, cross-correlating various arti- 
ficially broadened template spectra with DoAr 21 to find the template that gives the highest 
correlation peak, and we consistently find vsini ~ 80 km s _1 . However, our spectra cover 
only two windows of a few days each, so it is possible that both our 2002 and 2003 ob- 
servations were obtained at orbital phases with relatively large velocity differences. If we 
assume an equal- mass binary, the evolutionary tracks in Figure [T] give masses of ~ 1.8 M Q 
for each star, or a total sys tem mass of ~ 3.6 M . With an assumed semimajor axis of 
1-2 AU (ILoinard et al.ll2008l ). the relative orbital velocity of the two stars in a circular orbit 
would be ~40-60 km s _1 , which is of the right mag nitude to account for th e difference in 
line broadening observed in our spectra and those of IMassarotti et al.l (120051 ). The inferred 
orbital period of 0.5-1.5 yr from these orbital parameters is also consistent with our two 
observations (separated by ten months) falling either half or a whole orbital period apart, 
so that the same velocity difference would be observed. 



DoAr 21 has at times shown Ha in emission in the past (lHarolll949i ; iDolidze fc Arakelyan 



19591 ; lHidajatlll96ll ). while other observations have shown Ha absorption. Our observations, 
shown in Figure [TJ show a combination of both. There is a central photospheric absorption 
line, but our high spectral resolution reveals the presence of b road emission wings. I n 
addition, comparison with the Kl V spectral standard HD 13445 (IKeenan &: McNeill Il989l ). 
artificially broadened to vsini = 80 km s _1 , shows that DoAr 21's Ha absorption line 
is substantially filled in with emission. After correcting for the underlying photospheric 
absorption, the emission component we observe in DoAr 21 has an emission equivalent width 
of 1-2 A. 



iBarv et alJ (|2003f ) suggest that the inclination DoAr 21 is greater than 55°. If this is the case, then the 
equatorial rotational velocity is greater than 100 km s _1 . 



- 19 - 




-1000 -500 500 1000 

Velocity (km s _1 ) 



Fig. 7. — High-resolution optical spectra of DoAr 21 in the vicinity of the Ha line. Top: 
The dot-dashed line (red in the online edition) shows DoAr 21 on 12 April 2003; the solid 
line (blue in the online edition) shows DoAr 21 on 13 April 2003; the dotted line shows the 
Kl V standard star HD 13445, artificially broadened to vsini = 80 km s _1 . The Ha line in 
DoAr 21 clearly is filled in by a broad emission component that is present on both nights, 
but stronger on the second night. Bottom: The difference between the two DoAr 21 spectra, 
showing the broad emission component that appears the second night. The full-width at 
10% intensity is indicated. 
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Further, the line varies substantially from one night to the next. Figure [7b shows the 
difference between two nights, clearly showing the presence of a broad, variable emission 
component at Ha. The emission shown here is a lower limit to the true Ha emission line 
flux, since any emission component present on both nights is subtracted out. The FWHM 
of the emission is ~ 300 km s _1 , while the full width at 10% intensity is ~ 450 km s _1 . 



White fc Basril (120031 ) propose that stars with Ha full width at 10% intensity of > 270 km 
s _1 should be considered classical T Tauri stars, arguing that the line velocity width is a 
better accretion diagnostic than the more commonly used equivalent width; they observe a 
range of line widths from 270 to 670 km s _1 in accreting sources. By this measure, DoAr 21 
could be considered to be a classical (accreting) T Tauri star, and indeed iBary et al.l (120021 ) 
suggest that DoAr 21 is still accreting material from a disk based on their detection of H2. 

While it is possible that DoAr 21 could have some low-level accretion, its Ha emission 
is also consistent with a lack of accretion, since some star s that are no t thou ght to be 
accreting show Ha behavior similar to that seen in DoAr 21. iMontes et al.l (119971 ) find that 
some chromospherically active binary stars show Ha emission that is well-fit by a broad 
component and a n arrow component, with t he broad component having FWHM values from 
133 to 470 km s _1 . iFernandez et al.l (120041 ) find that the quiescent (non-flare) spectrum of 
the weak-lined pre-main-sequence triple system V410 Tau has Ha that is just filled in, not 
overtly in emission, while during a flare the Ha equivalent width is 27 A. After subtraction 
of a photosphere of the correct spectral type from the quiescent spectrum, they find that the 
residual Ha emission line (equivalent width ~ 1 A, similar to the emission component from 
DoAr 21) has a broad component with a FWHM of 300-400 km s _1 . During a flare, the 
H/? line is observed to be in emission with FWHM of 730 km s _1 , with its width declining 
as the flare decays. V410 Tau shows no evidence of infrared excess or other indicators 
of accretion, though definitive detection of a very weak infrared excess might be difficult 
given the shortage of resolved infrared measurements of the three stars (all within 0'.'3; 
White fc Ghe^ 200 l|) in the syst em. There is no det ectable emission in t he CO fundamental 
jNaiita. Carr. fc Mathieull2003h nor H 2 FUV lines faerczeg et aDl2006h . With no evidence 
of either a gas or a dust disk around V410 Tau, the Ha seems more likely to be related to 
the flares than to accretion. Thus, broad Balmer emission lines are not by themselves an 
unambiguous indicator of accretion. 

We revisit the question of what the Ha and H2 emission tell us about DoAr 21 's circum- 
stellar material in Section HI Here we simply note that the chromospheric activity implied 
by one interpretation of the Ha emission, and the need for a strong UV flux to excite the 
H 2 and PAH emission, both imply a need for a better understanding of the high-energy 
emission from DoAr 21. Thus, we now turn to an examination of our new high- resolution 
X-ray spectrum of DoAr 21. 
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3. X-ray emission 



DoAr 21 is the brightest X-ray source in the p Oph cloud core A, with an X-ray lu- 
minosity of nearly 10 32 ergs s _1 , a hard thermal spectrum, and large X-ray flares seen in 
several archival observations at a rate of nearly one per day. It s X-ray e mission has been 



obser ved at lo w resolution with Ei nstein ( IMontmerle et al. 



19951 1 . ASCA flKovama et al.lll994f ). and Chandra (llmanashi et alJbOOa iGagne et al.ll2004f ). 



1983jV ROSAT flCasanova et al. 



The new Chandra grating spectra we present here have vastly superior spectral reso- 
lution compared to any of these previous datasets, enabling us to extract temperature and 
abundance information with more reliability and also, for the first time, to examine density 
diagnostics. Given the high rate of flaring seen in previous observations, it is not surpris- 
ing that we detected a large flare in our 94 ksec Chandra observation. Thus, we can bring 
the same temperature and abundance diagnostic tools to bear on the pre-flare, flare, and 
post-flare portions of the new dataset. 

The High Energy Transmission Grating Spectrometer (HETGS) has two grating arrays: 
the Medium Energy Grating (MEG), with a FWHM resolution of 2.3 mA ; and the High 
Energy Grating (HEG), with a resolution of 1.2 mA (ICanizares et al.ll2005l ). Both grating 
arrays operate together, with the dispersed spectra (first, second, and third orders) as well 
as the zeroth order spectrum recorded on the ACIS CCD array. We used standard CIAO 
(v.3.3) tools to extract the dispersed spectra, and to create observation-specific spectral 
response matrices and effective area tabulations, as well as to create light curves and also to 
extract spectra for each of the three subsets of the observation. Of the dispersed spectra, 
only the first order spectra have a significant number of counts. The zeroth order spectrum 
is severely piled up, so we restrict the analysis we report on here to the first order MEG and 
HEG spectra. We performed most of the spectral analysis using xspec v. 12.3. 

The goal of our X-ray spectral and timing analysis is to characterize the properties of 
the hot plasma on this very magnetically active pre-main-sequence star, and to compare 
its properties to those measured in other pre-main-sequence stars with high-resolution X- 
ray spectra. Its extreme youth and lack of obvious signatures of accretion make for an 
interesting contrast between DoAr 21 and th e strongly accreting T T auri stars that have 
been well-studied with the Chandra gratings (IGiidel fc Telleschil 120071 ). Although DoAr 21 
is not a classical T Tauri star, it does — as we have detailed in the previous section — have 
some circumstellar material in its vicinity. The observed X-rays, and the unobservable far- 
and extreme-UV emission associated with the X-ray emitting plasma as it cools, could have 
an important effect on this circumstellar material, perhaps being the dominant source of 
excitation for the PAH and H 2 emission. Thus, characterizing the X-ray properties of the 
star is important for understanding the physical conditions in the circumstellar environment. 
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In Figure [8] we show the MEG and HEG spectra (negative and positive first orders 
co-added in both cases), with emission lines identified and labeled. Although the Chandra 
gratings and detector have significant response to wavelengths above 30 A, interstellar at- 
tenuation (due to photoelectric absorption, which has a cross section that goes roughly as 
A 3 ) makes the spectra above 12 A nearly devoid of counts. These spectra are dominated 
by a strong bremsstrahlung continuum, which is indicative of plasma with a dominant tem- 
perature well in excess of 20 million K (so that atoms are mostly fully stripped and their 
associated line emission is weak). The presence of emission from high ion stages — up to 
helium-like Fe xxv and hydrogen-like Fe xxvi — also indicates very high plasma temper- 
atures. We analyze the temperature distribution in the plasma in detail below by fitting 
thermal emission models to the entire spectrum. 



3.1. Spectral Modeling 



We fit the MEG and HEG first order dispersed spectra simultaneously (but not co- 
added) over the spectral ranges where each has a significant number of counts. For the 
MEG this was 2-12.5 A, and for the HEG it was 1.5-11.5 A. We fit a two-temperature op- 
tically thin thermal e mission model (the bapec implementation of the Astrophysical Plasma 
Emission Code apec, ISmith et al.l l200ll ) that accounts for bremsstrahlung and line emis- 
sion from a plasma in statistical equilibrium. This model has four free parameters: the 
plasma temperature, the abundances (expressed as a fraction of solar), the emission mea- 
sure (EM = 4nd 2 j n e nftdV, proportional to the normalization of the model), and the line 
broadening, which is an ad hoc turbulent velocity added in quadrature to the thermal ve- 
locity i _of_^achJine_jn_Jhe_jrmdel. We a lso include interstellar attenuation, with cross sections 
from iMorrison &: McCammonl (119831 ). We used the \ 2 statistic with Churazov weighting 



( Churazov et al. 19961 ) to assess goodness of fit and to place confidence limits on the derived 



model parameters. 

A single temperature model does not provide a good fit, though the low resolution 
ASCA and Chandra ACIS data are adequately fit by a single temperature thermal model 
( llmanashi et al.ll2002l ). We did find a good fit when we used a two temperature bapec model 
with interstellar absorption. The best-fit model has temperatures of roughly 12 and 48 mil- 
lion K (MK) , with approximately five times the emission measure in the hotter component as 
in the cooler component. The abundances are sub-solar, and no significant line broadening 
is found, with a 68% confidence limit of a tU rb = 50 km s -1 , about one-third of the peak 
spectral resolution. We find an interstellar column density of slightly more than 10 22 cm -2 , 
which is completely consistent with the extinction of Ay = 6.2 magnitudes, given the con- 
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Fig. 8. — The entire usable portions of the MEG (top) and HEG (bottom) co-added (negative 
and positive) first order spectra of DoAr 21. The bin sizes are 2.5 mA for the HEG and 5 
mA for the MEG. The shape of the continuum is dominated by the wavelength-dependent 
effective area of the telescope, gratings, and detector, and by interstellar absorption, rather 
than by the intrinsic emission. Vertical dashed lines represent the laboratory rest wavelengths 
of all detected emission lines. The annotations above the top panel indicate the element and 
ionization stage of each line, while those above the lower panel highlight the Lyman-a lines 
of the hydrogen-like isoelectronic sequence and the principal transition of the helium-like 
isoelectronic sequence. 
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version between extinction and hydrogen column density of IVuong et al.l (120031 ). The model 
has a flux of 1.48 x 10~ n ergs s -1 cm -2 on the range 0.3-10 keV. Correcting for interstellar 
attenuation, the unabsorbed X-ray flux is a little more than twice this, corresponding to an 
X-ray luminosity of Lx = 5.44 x 10 31 ergs s _1 . The best-fit model parameters and their 68% 
confidence limits are listed in the first row of Table [2j This fit is formally good, although the 
two temperatures certainly are an approximation to a continuous distribution of tempera- 
tures. Furthermore, the confidence limits are based only on statistical errors on the data, 
and in our experience, for datasets with many bins (the data we fit here have 12,196 bins), 
formal confidence limits are unrealistically tight. 



To test t he robustness of this fit, w e refit the data using the ISIS X-ray analysis package 
(ver. 1.4.9-38; iHouck fc Denicolall2000l ). We used the APEC model with interstellar absorp- 
tion, and, as we did with the fits in xspec, found that a two temperature model was required 
to achieve a good fit. The temperatures of the best-fit components were both higher by about 
20%, and the emission measure weighting was even more skewed toward the hot component 
in this ISIS fit compared to the xspec fit. The interstellar column density was about 10% 
lower than in the xspec fit. There were s ome differen ces in the method used for the fitting 
in ISIS, including the use of the C statistic (]Cashlll979l ) and a somewhat different wavelength 
range (1.8-18 A for both the MEG and HEG data). Perhaps these 10-20% discrepancies 
between the results from the two model fitting programs are a more realistic representation 
of the parameter uncertainties than are the formal, statistical confidence limits. We will use 
the apec model fitting in xspec as the standard throughout the rest of the paper. 

The best-fit model to the full observation, with parameters listed in the top row of 
Table [2j reproduces all portions of the spectrum quite well, with few systematic deviations. 
The best-fit model is shown in Figure M along with the higher signal-to- noise MEG data. In 
general, the continuum is well fit and the lines are adequately reproduced by this model. 
The same is true for the HEG data. 



3.2. Time Variability and Spectral Analysis of the Pre-flare, Flare, and 

Post-flare Data 

A large (factor of three) flare was seen in the middle of the Chandra observation, with 
a rapid rise of about six thousand seconds, followed by a steady decay during which the 
count rate drops by less than a factor of two in about 18 ks. After this, the count rate 
remains steady, but at a level that is nearly twice the quiescent, pre-flare level. The temporal 
behavior is shown in Figure [TOj where we also display light curves for the hard and soft bands 



Table 2. apec model fits 



(keV) 



EM X 
'10 53 cm- 3 ) 



kT 2 
(keV) 



EM 2 
10 53 cm- 3 ) 



Abundance 
(Solar units 3. 



N H {ISM) 
(10 22 cm" 2 ) 



erg s" 



10 31 — ,-1 



total 


1.00^ 


-.04 
-.03 


6.82+^ 


4.16 ± .11 


30.7+^ 


0.39 ±.02 


1 1O+-03 
i - iJ -.04 


5.44 


pre-flare 


0.95^ 


-.05 
-.06 






2i.8i; : ? 


0.26 ± .02 


1 1Q+-03 
i -- Ly -.05 


3.49 


flare 


1.58j 


-.08 
-.10 


i8.i±|;l 


7.81+ 2 ; 22 


40.5±|;| 


0.431$ 


1.18 ± .05 


9.36 


post-flare 


0.80^ 


-.09 
-.05 


n-i±?:l 


4 4Q+-24 
^•^ y -.50 




0.31±-« 




6.78 



a Solar abundances from lAnders fc Grevessd (119891 ) 
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separately. Significant hardening is seen during the flare, which dissipates through the flare 
and post-flare phase but never returns to the pre-flare level. 

Above we have presented fits to the total spectrum collected during the 94 ks observa- 
tion. We next repeat the two-temperature thermal fitting on the pre-flare (0 through 48 ks), 
flare (54 to 70 ks), and post-flare (70 to 94 ks) spectra separately. The results are reported 
in Table EJ and the temporal divisions are indicated by the vertical lines in the first panel of 
Figure [TUl The temperatures and emission measures roughly double during the flare portion 
of the observation, with somewhat larger relative changes in the hotter component of the 
two-temperature model. The flare plasma temperature is in excess of 90 MK. The post-flare 
plasma also has a temperature distribution that is elevated with respect to the pre-flare 
plasma, indicating that some reheating is occurring. The absorption column densities de- 
rived for each section of the observation are completely consistent with each other, which 
is consistent with an interstellar origin to the attenuation, as we would expect. We also 
note that the abundances show a marginally significant change among the three different 
sections of the observation, in the sense that the hotter sections — the flare and to a lesser 
extent the post-flare section — have higher abundances. Although this result may represent 
a real change in the coronal plasma abundances during the flare, it may also be affected 
by the simplified two-temperature assumption in the model. The abundance parameter in 
the emission model is largely controlled by the line-to-continuum ratio, and the strength of 
the metal emission lines relative to the bremsstrahlung continuum is also affected by details 
of the emission measure distribution at high temperatures, where metals are generally fully 
ionized and contribute little overall line emission. 

Another demonstration of the spectral changes associated with the flare is provided by 
the iron K-shell line emission, which is sensitive to the hottest plasma. The He-like Fe xxv 
line complex near 1.86 A has an emissivity that peaks near 60 MK, while the H-like Fe xxvi 
line at 1.75 A has an emissivity that peaks at temperatures somewhat above 100 MK. There 
is very little emission in the He-like feature and none detected in the H-like feature during the 
pre-flare section of the observation. But the iron lines get significantly brighter during the 
flare and maintain significant emission levels during the post-flare phase as well. Not only 
do the iron features increase in overall intensity during the flare, but the relative strength of 
the H-like feature compared to the He-like feature also evolves in the sense expected from 
the temperature changes seen in the apec model fitting to the entire spectrum. 
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Fig. 9. — The best-fit two-temperature thermal emission model is superimposed (in gray; 
red in the on-line edition) on the MEG data, with fit residuals shown below. 
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Fig. 10. — Left: A light curve with 1000-second bins formed from all counts in the dispersed, 
first-order spectra (both MEG and HEG). Right: Light curves made from all counts with 
wavelengths longer than 6.8 A (open circles) and shorter than 6.8 A (filled triangles). The 
hardening of the emission during and after the flare is evident. 
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3.3. Density-Sensitive Emission Line Ratios 

The emission lines in the spectrum are relatively weak, but we are able to analyze the 
strengths of two important helium-like line complexes and use them as plasma diagnostics. 
Specifically, the forbidden-to-intercombination line ratios of helium-like ions are sensitive 
to density, as collisions de-populate the metastable upper level 3 S of the forbidd en line 



and populate the upper level 3 P of the intercombination line (IGabriel fc Jordan! 1 19691 ). Thus 
higher densities decrease the forbidden-to-intercombination (f/i) ratio (with different critical 
densities for each element). This evidence for high densities is seen in some acc reting T Tauri 



stars , but not in naked T Tauri or magnetically active main-sequence stars (IKastner et al. 



20021 ; iTelleschi et al.l 120071 . and references therein) . 



We measured the line intensities by first fitting the continuum near each line complex, 
and once the continuum level was established, fitting the line complex itself with a three- 
Gaussian profile model on top of the best-fit continuum level. The adjustable parameters 
of the three-Gaussian model are an overall normalization and the ratios Q = (/ + i)/r and 
1Z = f/i, where r represents the flux in the resonance line. The f/i ratios for Si xm and 
S XV are 5.41 2 ,' 2 an d l-9*g 'f, respectively. None of the lower atomic number elements have 
enough signal-to-noise in their helium-like emission complexes for a measurement to be made. 

For the Si xm complex, the forbidden line is blended with the Lyman-/? line of Mg 
XII. Based on the measured Lyman-a flux of Mg XII and the relative emissivities of the 
two components of the Lyman series in the apec model, we can attribute up to 20% of 
the measured flux in this feature to the Mg line. This makes the ratio 1Z = f/i — A.^t^'t- 
The models predict f /% < 2.3-2.5, with the largest line-ratio values (the "low-density limit") 
expected for all electron densities less than or equal to a certain limiting valuejfl The 68% 
lower confidence limit on the f/i ratio from our data is 2.1, which limits the density to less 
than 4-9 x 10 12 cm -3 , depending on the model of lZ(n e ) that is employed. The MEG and 
HEG measurements of the line complex and the theoretical models from which we derive the 
upper limit on the density are shown in Figure [TT1 



The low-density limit refers to the limiting (high) value of the f/i ratio that holds when the density is 
low enough that collisional excitation out of the metastable excited state of the forbidden line is unimportant 
compared to spontaneous emission to the ground state. We give a range of value s for the low-density f/ i 
limit for each ele ment. These represent th e range of values found from PrismS pect (jMacFarlane et al.ll2004l ) 
calculations, from Blumenthal et al. ( 1972 ). and from Porquet fc Dubau ( 2000f) . If the actual electron density 
were much above the low-density limit, then the forbidden line would be weakened and the f/i ratio would 
be reduced. 
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Fig. 11. — The Si xm complex in the HEG (left) and MEG (center) shows a strong forbidden line (at 6.74 A) and a 
weak intercombination line (at 6.68 A), consistent with the low-density limit. The wavelengths of the three components 
(resonance, intercombination, forbidden from short to long wavelengths) are indicated by the dashed vertical lines, 
and the solid line shows the best simultaneous fit to the M EG and HEG data . Righ t: Models of the f H ratio of Si 



xm a s a function of electron den sity. The models are fr om iBlumenthal et al.l (I1972L solid curve), iPorquet fc Dubau 
(120001 . dotted), and PrismSpect (IMacFarlane et al.l 12004 dashed). The horizontal dash-dot line represents the 68% 
lower confidence limit of f/i = 2.10. Depending on which of the three models are used, this value corresponds to an 
electron density of n e = 4-9 x 10 12 cm -3 . 
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For S XV the forbidden-to-intercombination ratio we measure is consistent with the 
low-density limit for that complex of f/i = 2.0 to 2.1. The 68% lower confidence bound 
on f/i for S XV corresponds to a density of 10 14 cm -3 (which represents an upper limit 
to the density). That the helium-like complexes on DoAr 21 have strong forbidden lines, 
consistent with being in the low-density limit, is in marked contrast to altered f/i ratios 
in accreting classical T Tauri stars (CTTS). However, the altered f/i ratios in CTTS are 
generally seen in lower Z elements' He-like complexes, such as Ne IX and O VII, which 
have lower critical densities than Si xm. The lower Z features seen in those CTTS, which 
are at longer wavelengths, are not visible in the Chandra spectrum of DoAr 21 because of 
the large amount of interstellar attenuation. Nonetheless, the available density diagnostics 
for the X-ray emitting plasma of DoAr 21 are completely consistent with those seen for 
active, non-accreting coronal sources such as weak-lined T Tauri ( WTTS) and ZAMS stars. 
Specifically, the upper limit we find for n e from the Si xm complex is slightly lower than 
the densities in the accreting T Tau ri star TW Hya, where n e = 10 is found from the 
Ne IX complex (IKastner et al.ll2002l ). and our upper limit is also consistent with values ~ 2 
orders of magnitude lower than that, typical of coronal emission seen in RS CVn stars. 



4. Discussion 

Our X-ray observations show DoAr 21 to be a very active, hard, coronal source, with 
a large X-ray luminosity and flare rate. At the same time, our new infrared observations of 
DoAr 21 show the presence of extended, irregularly distributed circumstellar material. In 
light of this combination, we now consider the relationship between DoAr 21 's radiation and 
its circumstellar environment, and what it can tell us about DoAr 21's evolutionary state 
and more generally about disk evolution around young stars. 

Much of what we can learn from these data rests on where exactly the circumstellar 
material is located, and how it is related to DoAr 21. Is it the remnant of a disk that is in the 
final stages of being cleared out? Or is it nearby but not necessarily gravitationally-bound 
material that is illuminated by DoAr 21 as it emerges from or passes through the cloud? 

We divide the discussion into two major parts. First, we consider the emission mecha- 
nism and energetics of the observed infrared emission (including the H 2 and PAH emission) 
from the circumstellar material; as a corollary to this we evaluate further literature data 
bearing on DoAr 21's level of stellar activity, and we use it to estimate the far-ultraviolet 
(FUV) radiation field. We then discuss possible origins for this circumstellar material, fo- 
cusing in particular on whether or not it is a disk, i.e. gravitationally bound, centrifugally 
supported material. 
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4.1. Energetics and excitation of the circumstellar material 

The most surprising result from the new observations presented here is the large spatial 
extent of the mid-infrared emission from DoAr 21, which is puzzling to find around a cool 
star, and one with a very low millimeter flux. What this discovery has in common with the 
previous surprises about DoAr 21 — the detections of PAH and H 2 emission — is that all of 
these observations are more commonly associated with more luminous, hotter stars. Here 
we explore whether it is possible to construct a model that explains all of the data, and we 
focus in particular on trying to understand the influence of the UV and X-ray emission from 
DoAr 21 on its circumstellar material. Because FUV radiation is important for radiative 
excitation of both H 2 and PAH emission, we first discuss the FUV radiation field near DoAr 
21. 



4-1.1. The FUV radiation field near DoAr 21 



The ultraviolet flux from DoAr 21 is poorly known. It was not detected by IUE 
( TValenti et al.l 120031 ); examination of the archived IUE spectrum places a rough upper limit 
of ~ 10~ 15 erg cm -2 s" 1 at A = 2400-3200 A. Given the large extinction to DoAr 21 (roughly 
13 magnitudes at A = 2400 given our adopted visual extinction), however, the star could 
have an ultraviolet flux that is several orders of magnitude above its photospheric flux and 
it would not have been detected by IUE. Given this situation, we take two approaches to 
estimating the level of UV radiation in the vicinity of DoAr 21. First we consider the UV 
background in the region due to the nearby B2 star HD 147889. Then we use the available 
data on DoAr 21 's chromospheric and coronal activity to estimate the UV emission from its 
transition region. 

The B2 V star HD 147889, at a projected distance of 0.4 pc from DoAr 21, is responsible 



for ex citing PAH and H 2 emission from much of the surrounding nebulosity (IHabart et al. 



20031 ). To determine the UV flux from HD 147889 arriving at the position of DoAr 21 and 
its surroundings, we integrated the flux from a synt hetic spectrum ba s ed on a spherical, 
hydrostatic, class (ii) NLTE model atmosphere from lAufdenberg et al.l (119991 ) with T e g- = 
22000 K and \ogg = 3. 9, assuming L = 5 300 L & for HD 147889 and a distance of 0.4 pc 
between the two stars (ILiseau et al.lll999l ). At this distance, the ultraviolet flux reaching 
DoAr 21 from HD 147889 (assuming no attenuation) is roughly 0.75 erg cm -2 s _1 in the 
bandpass 6.0-13.6 eV (A = 912-2066 A). This corresponds to Go ~ 470, where Go is the flux 
in this wavelength interval expressed in units of a mean interstellar radiation field of 1.6 x 10 -3 
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erg cm " 2 s" 1 JHabinej [l968lR. In contrast, the FUV field from DoAr 21's photosphere 
(estimated using a Kurucz model atmosphere and the stellar parameters given above) is 
Go ~ 220 at a distance of 10 AU from the star and falls with distance from DoAr 21. As 
noted above, however, given the high level of activity from DoAr 21, there is likely to be 
significant UV emission from the transition region. 

To estimate the FUV flux from DoAr 21, we use as a template HD 283572 (V987 Tau), 
a weak-lined T Tauri star with a similar mass, effective temperature, and X-ray properties. 
HD 283572 is a single G5 star in the Taurus- Auriga star-formin g region, and it has an X-ray 



luminosity of Lx = 1-1 x 10 31 erg s 1 in the 0.3-10 keV b and (Telleschi et a 



2007, scaled 



to the 129 pc distance measured from VLBA parallax by iTorres et al.l 120071 ). DoAr 21's 
luminosity in this energy band is Lx = 3.5-9.4 x 10 31 erg s" 1 (Table [2]), taking the range 
from pre-flare to peak flare values. HD 283572 is near DoAr 21 on the HR diagram, although 
with a slightly hotter photospheric temperature and older inferred age; it is the second blue 



square from the left in Figure [Q It has no eviden ce for circumstellar dust (IFurlan et al.ll2006 



Kundurthv et al.l 120061) or IT (iBary et al.l 120031 ) . It shows a hard, flaring X-ray spectrum 
( Favata. Micela. fc Sciortino 19981 ); fit ting of the spectrum yields a temperature distribution 



similar to what we find for DoAr 21 ( 



Telleschi et a. 



radio emission flPhillips. Lonsdale, fc Feigelson!ll99ll ). 



20071 ). Like DoAr 21, it shows strong 



These similarities make HD 283572 a good template to co mpare with DoAr 2 1, along 
with the crucial difference: it has a relatively low Ay = 0.4 mag (IFurlan et al.l 120061 ) . Thus, 
the FUV emission is not attenuated nearly as much as for DoAr 21, and the star was clearly 
detected b y IUE. We integrated the de-reddened FUV flux in the co-added HD 283572 IUE 
spectrum (jValenti. Johns-Krull. fc Linskyl 120001 ) in the range A = 1150-2066 A. The Ly 
a line is clearly contaminated by geocoronal emission, and in any case would be strongly 
absorbed by interstellar H I; we rejected the flux from A = 1200-1232 A. Since the A = 912- 



1150 A flux was not observed by IUE, we used the results of Alexander. Clarke, fc Pringld 



(12005), who calcula t ed UV spectra using detailed differential emission measure distributions 



from iBrooks et al.l (120011 ) for several T Tauri stars and found that the A = 912-1150 A 
flux was roughly 1/6 of the A = 912-2000 A flux. Thus, we scaled the de-reddened IUE 
1150-2066 A flux by 1.2 to estimate the flux in the full 912-2066 A interval. We used the 
measured distances of HD 283572 and DoAr 21 (both with high-quality VLBA parallaxes) 
to determine the scaling between flux and luminosity for both stars. Finally, we scaled the 
HD 283572 FUV luminosity by the two stars' relative Lx values to estimate the DoAr 21 
FUV luminosity, assuming that the FUV luminosity would be directly proportional to the 



6 This i s the same quantity referred to as Gfuv bv lNomura k Millar! (|2005h and similar to the quantity 
X used bv lHabart eTal] |2003h 
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X-ray luminosity. 

Based on this analysis, we find Lfuv ~ 2.8-7.4 x 10 31 erg s _1 . However, this analysis 
does not account for flux in the Ly a line, which can be very strong in stars with active 
chromospheres. We estimated the Ly a flux usin g the scaling relations between Ly a and 



X-ray surface flux derived by [Wood et al.l (120051 ). who measured Ly a from a sample of 



nearby late-type stars at high spectral resolution with STIS on HST, and thus were able 
to remove geocoronal emission and correct for interstellar absorption. The range of X-ray 
fluxes from DoAr 21 gives a range of Ly a surface fluxes of 0.8-1.3 x 10 7 erg cm -2 s^ 1 , and 
thus a total Lpuv ~ 3.6-8.7 x 10 3 1 erg s " 1 . We caution, however, that there is only one star 



(Speedy Mic) in the IWood et al.l (120051 ) sample with an X-ray surface flux as large as that 



from DoAr 21, and so the scaling relations are not well-constrained in this regime. 

Under these assumptions, we estimate that the G value produced by DoAr 21 's FUV 
flux is in the range of 470-1100 at l'/l (134 AU) from the star, the distance of the bright ridge 
in our images. Thus, the estimated FUV flux from the star is equal to or somewhat greater 
than the external FUV flux at this distance. The total FUV flux (from the combination of 
DoAr 21 's transition region and the B2 star background) at the position of the bright ridge 
is then Go ~ 940-1600; the internal and external FUV fluxes become equal at ~ 130-210 
AU from DoAr 21. If there is any attenuation of the external FUV flux in the > 0.4 pc 
between HD 147889 and DoAr 21, then DoAr 21's contribution to the total flux is even more 
dominant. 

Having estimated the FUV flux near DoAr 21, we can now examine the effect this 
flux would have on nearby material. In the sections below, we will examine in turn the 
continuum (dust) emission, the PAH emission, and the H2 emission. Consistent with the 
FUV flux estimates above, we will show that the emission shows signatures of being excited 
from the inside out, i.e. by DoAr 21 rather than externally by HD 147889. Before turning 
to a discussion of the details of each type of emission, we close this section with a simple 
check of the energy balance in the circumstellar environment. We measured the flux from 
the circumstellar material by taking the fluxes in the largest aperture listed in Table [1] and 
subtracting the fluxes in the smallest aperture to remove the stellar contribution. Integrating 
from 10 fim to 20 /im, we find the luminosity of this material to be roughly 9 x 10 31 erg 
s -1 , the same order of magnitude as our estimate of Lfuv- Given the uncertainties in our 
flux estimate and the complexities of the various heating, cooling, and emission processes 
in the environment we do not suggest that there should be precise agreement here, but the 
order-of-magnitude agreement does indicate that DoAr 21 is a plausible source for powering 
the observed emission, which we now examine in more detail. 
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4-1.2. Dust emission 

In order to study the continuum emission excess from DoAr 21, which we here assume 
arises from thermal emission from dust grains, we need wavel e ngth regions without signifi- 



cant spectral features. The spectrum taken by lHanner et al.l (119951 ) shows that the region 
encompassed by the 10.4 jum filter from our Gemini observations is free of emission or ab- 
sorption features; in particular, there is no evidence of a silicate feature. While we do not 
have a corresponding spectrum in the 18.3 /im region, this part of the spec trum generally ha s 



not shown significant features in Spitzer IRS spectra of other young stars (IGeers et al.ll2006l ). 



Some galaxies and star-forming regions with PAH e mission do show features at 17.8 /im and 



18.9 /xm, albeit on top of a significant continuum (jSellgren. Uchida. fc Werner! 120071 ). We 
take this spectral region to sample primarily continuum emission; if there is a contribution 
from some PAH emission in this bandpass, the effect would be to raise the dust temperatures 
derived below. 

Using our Gemini images, we integrated the flux in these two bandpasses in a region 
that is off the stellar position and thus probes only circumstellar flux. This yielded fluxes of 
197 and 603 mJy at 10.4 and 18.3 /im, respectively. A single-temperature blackbody fit to 
these two fluxes yields a temperature estimate of roughly 215 K for the emitting material. 
While this estimate is obviously an oversimplification of the actual situation, the general 
result of a steeply-rising flux from 10-18 /im is robust, indicating that the emitting material 
has a temperature in the hundreds of degrees K. For comparison, blackbody grains that 
have reached thermal equilibrium with the star's radiation field, given DoAr 21 's luminosity 
of 11.7 L Q , would have a temperature of 44 K at 135 AU (l'.'l) from the star and 31 K at 
270 AU from the star. At these temperatures, the 18.3 /im to 10.4 /im flux ratio would be 
10 5 -10 7 , rather than the observed value of ~ 3. 

While the derived dust temperature is far out of equilibrium, it is quite consistent with 



the e xpected behavior for very small grains (VSGs) excited by UV photons (IPuget fc Leger 



19891 ). Grains with radii of a few nm can be transiently heated to temperatures of several 
hundred K by single UV photons, and an ensemble of such grains can contribute signifi- 
cant infrared e mission. Given that such small grains may form a continuum of sizes with 



PAH molecules ( iTielensI 120081 ) . it is perhaps not surprising to find them around DoAr 21. 
Nevertheless, it provides additional evidence for the importance of the FUV radiation in ex- 
citing the circumstellar material, and it is another respect in which DoAr 21's environment 
is different from a typical circumstellar disk. 
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4-1.3. PAH emission 



We have presented our new images in the 8.6 and 11.3 /im PAH filters above; here we 
briefly recap existing published and archival PAH data for DoAr 21. 



Figure [6] shows the mid-infrared spectrum of DoAr 21 from lHanner et al.l (119951 ) along 
with the mod el photospher e. Th e data show a clear emission feature at A = 11.3 /zm, 
attributed by lHanner et al.l (119951 ) to emission from PAHs. The data also show a rising 
slope shortward from A = 8.5 to where the spectrum ends at 7.9 /im, consistent with the 
presence of a 7.7 /im PAH emission feature. 

DoAr 21 was also observed by ISO with the circularly-variable filter (CVF) on ISOCAM, 
yielding a A/AA ~ 50-80 spectrum from 5.1- 16.6 ittm of DoA r 21 a nd its surroundings 
with 6" pixels. These data were reprocessed by iBoulanger et al.l (120051 ) as part of the final 
reprocessing of all ISOCAM CVF data, and we retrieved the reprocessed data from the ISO 
archive. The data show a strong point source at the position of DoAr 21, as well as a bright 
ridge of emission located 30"-90" NE of DoAr 21 and runni ng NW-SE; this ridg e is clearly 
visible in the Spitzer 24-/im image (Fig. [3j see also Fig. 3 of iPadgett et al.ll2008l for a larger 
view) and can also be seen faintly in the 2MASS K-band images. All pixels in the ISOCAM 
3f3 x 3'.3 field of view show a rich PAH spectrum, including emission features at 6.2, 7.7, 8.6, 
11.3, and 12.7 jum. 

Unfortunately, it is not possible to determine a reliable spectrum for DoAr 21 itself from 
these data. Due to internal reflections in ISOCAM when used with the CVF, the spectrum of 
DoAr 21 is contaminated b y "ghosts" , reflections of stray light from elsewhere in the field of 
view ( lOkumura et al.lll998l ). This is revealed by sharp discontinuities in the spectrum around 
9 /im, the transition wavelength between CVF1 and CVF2. The ghosts are significantly 
worse in CVF2, and this discontinuity (which is clearly present in a spectrum extracted a t 
the position of DoAr 21) is a signature of stray light contamination (IBoulanger et al.l 120051 ). 
We suspect that the contamination arises from the bright ridge of emission to the NE of 
DoAr 21, which is of the same order of magnitude in brightness at these wavelengths as 
DoAr 21 itself. Because the background emission is diffuse and spatially variable, and the 
ghost size and strength are functions of wavelength, it is not possible to reconstruct the 
uncontaminated spectrum of DoAr 21 from these data. 

One of the most striking features of our new PAH-filter images of DoAr 21 is the spatial 
variation of the 11.3 fim emission, with similar variations seen at 18.3 /im. The emission does 
not decrease smoothly with distance in the region near the star, but rather it is significantly 
brighter in a knot and partial arc beginning roughly l'/l to the north and west of the star. 
In contrast, the 8.6 fivcv emission is much smoother, with about the same brightness in the 



-36- 



region near the star as in the region that appears as a bright arc in the 11.3 /im image. 
Thus, it appears that the [8.6]/ [1 1.3] ratio may be higher near the star. However, the 8.6 
(Ltm emission is much fainter overall than the 11.3 /im emission, and the star is brighter, so it 
may be that any variation in the fainter 8.6 and 10.4 /im emission blends in with the wings 
of the stellar PSF. 

If the apparent change in the [8.6]/[11.3] ratio is real, it is a useful diagnostic of the 
ionization state of PAH molecules; in laboratory spectra of PAH molecules, the 11.3 /im 
CH out-of-plane bending mode weakens consider ably in ionized PAHs while the 8.6 £tm CH 



in-plane bending mode strengthens somewhat (lAUamandola. Hudgins. fc Sandfordl Il999l ). 



In astrophysical sources, a decrease in the [8.6]/ ]11.3] ratio with d istance from a central 



i onizing source has b een observed in NGC 1333 (IJoblin et al.lll996l ) and in the Orion Bar 



f lGalliano et al.l 120081 ) 



Applying this interpretation to DoAr 21 leads to the conclusion that a significant fraction 
of the PAH molecules within about 1" of the star are ionized, suppressing the [8.6]/[11.3] 
ratio in this region. This would indicate that DoAr 21, not HD 147889, is the dominant 
source of FUV flux in its circumstellar environment, consistent with our estimates above of 
the relativ e FUV fluxes (Sec. 14. 1. ID . The first ionization potential for typical PAH molecules 



is 6-8 eV (jWeingartner fc Draine 



2001bl ). and thus there must be significant flux at energies 



> 6 eV near the star. 

As noted above, it may be that the radial variation seen in the bright 11.3 /im emission 
is simply difficult to detect in the fainter 8.6 /im emission due to the low surface brightness. 
One piece of evidence in favor of this interpretation comes from comparison of the different 
continuum wavelength bands. Since both the 10.4 and 18.3 jum bands are dominated by 
continuum emission from very small grains, they should show the same spatial variation. 
However, the bright 18.3 /im emission shows the same radial variation as the bright 11.3 /im 
emission, while the fainter 10.4 /im emission apparently does not, which might suggest that 
it is not seen simply due to low surface brightness, and that a similar effect might be present 
at 8.6 /im. If it is the case that the same radial variation is present in all bands, then it 
would represent a real change in surface density of material, with less material present near 
DoAr 21, perhaps due to photoevaporation. In either case, the azimuthal distribution of 
material around DoAr 21 is clearly asymmetric (Fig. [5]). 

While there is a sharp increase in the 11.3 /im emission at roughly l'.'l from the star, 
the brightness of this emission declines with radial distance beyond that point. The decline 
is roughly consistent with a falloff proportional r -2 , which is what would be expected from 
a uniform PAH and small grain surface density with declining excitation due to geometric 
dilution of the exciting radiation. Thus, the true extent of the circumstellar material may 



be much greater than that shown in our images, as already suggested by the HST and 
near-infrared images (Sec. I2.3T21) . 



4- 1-4- H2 emission 



While we did not make any observations of H 2 emission, the analysis of the VSG and 
PAH emission in the preceding sections bears directly on the interpretation of previous 
H2 observations, which we discuss here. Given the surprisingly large spatial extent of the 
emission seen in our images, we pay particular attention to what constraints existing data 
put on the spatial distribution of the H 2 emission. 



As noted above, iBary et al.l (120021 . 120031 ) detected emission in the H 2 1-0 S(l) line at 
A = 2.12 /im from DoAr 21. Their slit width was 0'.'8 and the seeing was approximately I'M 
(IBary et al.ll2002l ). In their nod-differenced images, they found no evidence of any extension 
of the emission along the E-W slit beyond the l'/4 seeing disk. From the resolved velocity 
width of the line, they also argued that the gas must reside in a rotating circumstellar disk. 



Bitner et al.l (120081 ) observed DoAr 21 in the 0-0 S(2) pure rotational transition of H 2 
at A = 12.279 /im. They detected emission in one observation using a l'/4 N-S slit with l'/4 
seeing on the IRTF, but they did not detect it in a subsequent observation with a (X'54 E-W 
slit with 0'.'5 seeing on Gemini 



The Bitner et al. ( 20081 ) observations suggest that the S(2) H 2 emission is not uniformly 
distributed around the star and/or is time- variable. If we attribute the difference between 
the two observations solely to the spatial distribution of material, the data imply that the 
emitting gas lies primarily at an angular separation greater than C//25 (30 AU), and that 
it may be more concentrated to the N and/or S of the star. iBitner et al.l (120081 ) detected 
larger fluxes with the wider IRTF slit than with the narrower Gemini slit for four out of 
the five sources observed with both telescopes, suggesting that spatially-extended emission 
for the pure rotational transitions is not uncommon. We note, however, that a slit at any 
orientation will intercept flux from a disk at all disk radii (though not all azimuths), so 
azimuthally-symmetric emission from gas in a disk, even at large radii, should contribute 
some flux even in the Gemini slit. We discuss this point further in Sec. I4.2[ where we propose 
a possible geometry for the emitting material that is consistent with all the data. 



7 FWHM seeing profiles and slit orientations are from M. Bitner, personal communication. The stated 
slit orientations are the direction of the long dimension of the slit, i.e. the direction in which any extended 
emission would be detected. 
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The suggestion that the H2 S(2) emission may lie preferentially to the N or S, combined 
with the emission geometry seen in our PAH images, suggests that the H 2 and PAH emission 
could be coincident, and other observations suggest that this is not an uncommon occurrence. 
This is perhaps unsurprising, given that both molecules are excited by UV photons. In the 
Orion Bar PDR, the H2 1-0 S(l) emission and 3.3 /im PAH emission are similarly distr i buted 
alon g the edge of the bar facing the ionizing O stars^] Closer to DoAr 21, Gomez et al. ( 2003 ) 
and lHabart et aJj (120031 ) conducted ground-based narrow-band imaging, showing that the 
SW edge of the bright bar lying NE of DoAr 21 ( seen in both the Spi tzer and ISO images) 
is also bright in the H2 1-0 S(l) line. Further, lHabart et al.l (120031 ) show that the PAH 
emission seen by ISO CAM at A = 5-8.5 /im in this region correlates very well in strength 
and position (albeit with modest spatial resolution) with the H2 1-0 S(l) emission seen in 
the ground-based images, and with the H 2 0-0 S(3) emission seen by the ISO SWS. 



An alternate interpretation of the lBitner et al.l (120081 ) observations is that the emission is 
time-variable and was below the detection threshold at the time of the Gemini observation. 
The 3<j upper limit in the Gemini observation is 1.25 x 10~ 15 erg cm -2 s -1 , a factor of 
2.6 below the flux level of the IRTF detection^ The X-ray flux level from DoAr 21 is 
seen to vary by a factor of three on short timescales (Fig. [TUT) ; if this translates into a 
corresponding variation in the FUV exciting the H 2 emission, then that emission could vary 
by similar factors. However, if we attribute the differ ence between the tw o observations solely 
to variability, this hypothesis would require that the lBitner et al.l (120081 ) detection happened 
to coincide with the peak of a flare, and that the H 2 flux in the Gemini observation is just 
below t he detection threshold . DoAr 21 has a relatively high X-ray flare rate of about one 
per day (llmanashi et al.ll2002l ). but the fraction of time the X-ray flux is more than a factor of 
2.6 above the quiescent level is small, probably less than an hour per day based on available 
data. 



4-1.5. Photometric and Polarimetric Variability 

Multiple measurements exist at wavelengths from 0.39 /im (U) through 4.8 /im (M). 
Variability is seen at all wavelengths. The range at V band in 15 measurements is 0.33 mag- 



3 While to our knowledge there is not a similar map of the 11.3 /im line in the Orion Bar, the 3.3 /im PAH 



line i s strongest in neutral PAHs and in general is strongly correlated with the 11.3 /im line (|Honv et al 
200lh . 



9 This limit is determined using the line width of 5.6 km s 1 measured in the IR TF observations (M . 
Bitner, personal communication), and thus it is slightly lower than the limit given in iBitner et al. (2008) 
based on an assumed line width of 10 km s _1 . 



-39- 



nitudes, with variati ons on timescales of less than one day, but with no apparent periodicity 
( IBouvier et al.lll988l ). Substanti al variability is seen in the ultraviolet; in two observations 
separated by less than an hour, IBouvier et al.l (119881 ) observed a brightening of 0.97 mag- 
nitudes in U, while B brightened by 0.03 magnitudes and V faded by 0.02 magnitudes. 
Two observations separated by 3 hours on the following night showed a similar variation, 
with U brightening by 0.42 magnitudes while B and V faded by 0.04 and 0.05 magnitudes, 
respectively0 Such a large increase in the UV emission without a concurrent increase in 
V cannot be due to rotational mo dulation (either of hot or cool spots), since that would 



increase the flux in all bands (e.g.. iRydgren Sz Vrballl983l ). These [/-band brightenings are 



simila r to flares seen on acti ve main-sequence stars in that they are much brighter at U than 
at V ( [Fernandez et al.ll2004l ). We note that both of these brightening events occur when the 
star is at the faint end of the observed U range, i.e. they are brightenings toward the mean 
value rather than above it. This may simply indicat e that the star is fla ring much of the 
time and that the observed mean value in the data of IBouvier et al.l (119881 ) is characteristic 
of the flares rather than the photosphere; this is consistent with the fit to the SED (Fig. [2]). 
The U and B magnitudes are positively correlated with each other (with a larger amplitude 
at U), but neither the U or B magnitudes nor the U — B color is correlated with V. 



The bluest observed U — B color for DoAr 21 in the photometry of IBouvier et al.l (119881 ) 
is U — B — 0.5, observed on two consecutive nights. Using our adopted exti nction, the de- 
redden ed color is (U — B)o ~ —0.6. Interestingly, the only other stars in the IBouvier et al. 
( 119881 ) sample that show U — B excursions to such blue values are classical T Tauri stars, 
though the availability of well-sampled time series U and B photometry for some of the 
most active WTTS (e.g., HD 283572) is limited. On DoAr 21, where we see no evidence of 
accretion, the [/-band excess probably signals Balmer-line and continuum emission result- 
ing from extreme chromosph eric activity and flaring, as seen on late-type flare stars (e.g., 
Lacv. Moffett. fc Evans! 1 19761 ) . Flares from such stars have a mean U — B of —0.9 ± 0.3 



( lMoffettlll974l ). consistent with the bluest color seen from DoAr 21. 



DoAr 21's f^-band emission is polarized at p ~ 2% (IMartin et al.lll992l ). The large 
visual extinction toward DoAr 21 suggests that this could be attributed to interstellar po- 
larization from overlying cloud material not closely associated w ith the star. Howe ver, the 
polarization is variable both in amplitude and in position angle (I Jensen et al.ll2004l ). Since 



1( lBouvier et al.l (|1988I ) do not give error bars on individual measurements. They quote typical errors of 
0.09, 0.06, and 0.04 mags for U, B, and V for a V — 12 star; DoAr 21 is V w 14. Observations of other 
stars with similar V magnitudes observed during the same nights as the U flares on DoAr 21 do not show U 
magnitudes that are substantially different from the mean values for those stars, suggesting that the observed 
changes at U in DoAr 21 are real, and not photometric errors. 
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a cloud far in the foreground would not be expected to change on short timescales, this vari- 
ability suggests a circumstellar origin for DoAr 21 's polarization; this is also quite consistent 
with the observed asymmetry of the circumstellar material, as scattering off an asymmetric 
distribution will give a net polarization in the unresolved light. Given the discovery reported 
here that the circumstellar material around DoAr 21 is quite extended (Sec. 12.3.21) . the ob- 
served variation in polarization percentage could be due to measuring with different aperture 
sizes that in clude different amou nts of the extended emission: the percentage polarization 
me asured by iJensen et al.l (120041 ) in a 4'.'3-diameter aperture is smaller than that measured 
by iMartin et al. Jl992l ) in a 7'.'8 aperture. 



4.2. A Disk or Not? 



What is strikingly clear from the observations presented above is that the high-energy 
emission from DoAr 21 has strong interactions with material in its circumstellar environment, 
even on scales of hundreds of AU. What is less clear, however, is the relationship of the 
observed material to DoAr 21, and specifically to what extent the extended emission is 
related to any orbiting, bound circumstellar (or circumbinary) disk that DoAr 21 has or 
once had. 

The definitive answer to this will await further, higher-resolution observations, in par- 
ticular spatially- and velocity-resolved mapping of the H 2 emission. Here, we close the 
discussion with an examination of evidence for and against a Keplerian disk around DoAr 
21. 



One of the primary arguments that iBary et al.l (120021 ) made for the H2 emission arising 
from a disk is that extended, non-disk material would have a narrower line width than the 
9 km s _1 FWHM width they observed for DoAr 21's 1-0 S(l) emission. While we examine 
this line width more quantitatively below, we first note that DoAr 21 has the narrowest H 2 



line o f any of the detections of h 



SCI) faarv et~alll2003l: Etohelaj 



from T Tauri stars and Herbig Ae/Be stars at both 1-0 



2003: Takami et al 



2004 



Barv et alil20o5 beck et~alll2008h and 0-0 S(2) faitner et al.ll2008h . 



Ramsay Howat fc Greaves! 120071 : 



To interpret the line width, we must take into account the fact that the observed line 
width is the con volution of the intrinsic line width with the instrumental profile. The 
Bary et al.l (120021 ) observations had a FWHM spectral resolution of 5 km s _1 . Assuming 
Gaussian profiles for both the spectral response of the instrument and the observed line, the 
intrinsic line width is then 7.5 km s _1 . If i nterpreted s o lely a s a thermal line width, this 
would imply a kinetic temperature of 2450 K. IBary et al. J2002h do not quote an uncertainty 
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for the measured line width; assuming a 1 km s _1 total uncertainty for the deconvolved line 
width (reflecting some uncertainty in both the measured line width and the spectral resolu- 
tion), the range of derived kinetic temperatures is then 1700-3300 K. The lower end of this 
temperature range is q uite consistent wit h temperatures generally found for 1-0 S(l) emit- 
ting gas. F or example , Barv et al.l (120081 ) quote a standard thermal excitation temperature 



of 2000 K. iBeck et al.l (120081 ) observed extended H2 around several T Tauri stars, and from 



the ratios of several lines they de rive LTE temperatures of 1800-2300 K for H 2 around T 
Tauri stars. iHerczeg et al.l (120041 ) model the fluorescent H 2 emission from TW Hya's disk 



and find kinetic temperatures of 2500^5qq K for the H2. While the gas in some of those 
sources may be shock- heated, which is unlikely for DoAr 21, these observations show that 
H2 can survive and emit strongly at the tem peratures necessary to ex plai n the observed 



l ine w idth in DoAr 21. Similarly, the models of iNomura &: Millar! (120051 ) and iNomura et al. 



(120071 ) predict that the bulk of the H 2 1-0 S(l) emission from disks irradiated by X-rays 
and UV comes from gas with temper atures of 1000-2000 K, and similar results are found 
in models of P DRs (lAllers et al.ll2005r). Collisiona l dissociation of H 2 at these temperatures 
is insignificant (IMartin. Schwarz. fc Mandyl 119961 1 . Thus, the H 2 temperature necessary to 
explain DoAr 21's observed line width solely with thermal broadening is consistent both with 
theoretical expectations and with observations of other sources. 

Theoretical models show that it is plausible that DoAr 21 can heat gas to these tem- 
peratures at large distances from the star. The FUV flux that we estimate for DoAr 21 
combined with HD 147889 (Sec. KTJh corresponds to G « 940-1600 at a distance of 130 
AU (l'.'l) from DoAr 21, the location o f the b right PAH emission ridge. This is similar to 
the FUV flux in the INomura k Millar! (j2005h models at radii of ~ 40-100 AU, and thus 
we might expect conditions at 130 AU from DoAr 21 to be similar to those at 40-100 AU 
in the models. This range of radii extends beyond the ~ 30 AU distance at which the H 2 
1-0 S(l) emission peaks in those models, but it is also the case that DoAr 21 has an X-ray 
luminosity that is 50 times larger and an X-ray spectrum that is much harder than that as- 
sumed in these models (which base their assumed X-ray emission on that of TW Hya), thus 
contributing an additional source of heating. Thus, while the models do not apply directly 
to DoAr 21, plausible scaling from them suggests that H 2 at > 100 AU from DoAr 21 could 
have T > 1000 K, necessary to produce the observed line width with thermal motions, and 
could emit at the observed level, since the predicted H 2 flux in both the 1-0 S(l) and pure 
rotational lines in these models is similar to that observed for DoAr 21. 



To summarize this part of the discussion, the primary argument of iBary et al.l (120021 ) 
for the presence of a disk around DoAr 21 is essentially that the only way to produce the 
observed H 2 1-0 S(l) line width is from rotation in a disk. However, we have shown that the 
observed width is a plausible thermal line width for H 2 illuminated with a strong FUV and 
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X-ray flux, and that the flux from DoAr 21 is sufficient to produce the necessary heating at 
distances of > 100 AU from the star, consistent with the position of the observed PAH and 
VSG emission. While the match between this expected thermal broadening and the observed 
line width does not rule out the gas also having some bulk motion, for example from rotation 
in a disk, it does show that rotation is not required to explain the data. Indeed, the observed 
line width sets some significant limits on the w rbit sin i of the gas. The expected Keplerian 
velocity of gas orbiting a 2.2 M & star at 130 AU is 3.8 km s -1 . Given the uncertainties in the 
range of possible temperatures for the emitting gas, such an orbital velocity is likely allowed 
by the data for any sini. Orbiting gas at significantly smaller radii, however, would require 
values of sini significantly less than one. 



As noted above, the two measurements of H 2 from DoAr 21 by lBitner et al.l (120081 ) gave 
different results, suggesting that the bulk of the 0-0 S(2) emitting gas lies to the north or 



south of the star. However, the 5.6 km s 1 FWHM o 



not coincident with the gas detected by iBarv et a. 



the S(2) line traces cooler gas, so the iBitner et aL 



that is somewhat farther from DoAr 21 than the iBary et al.l ((20021) observations 



: the S(2) line suggests that this gas is 



(120021) in the 1-0 S(l) line. In general 



(200 81) observat ions may be tracing gas 



We have presented above a large body of data on DoAr 21, some of it seemingly con- 
tradictory. To conclude this discussion, we present here a broad picture of DoAr 21 's cir- 
cumstellar environment. This is not the only possible scenario, but we suggest that it is 
consistent with all of the available data, and it is testable with future observations. In this 
scenario, DoAr 21 no longer possesses any substantial orbiting circumstellar material, either 
in the form of gas or dust. Either the system never had much of a disk (possibly as the 
result of the formation of the binary companion) , or the extreme FU V and X-ray luminosity 
have photoevaporated the disk. The excess infrared emission, as well as the H 2 emission, 
comes from the extended material seen in our Gemini images, at hundreds of AU from the 
star. This material may be molecular cloud material that is not directly associated with the 
formation of DoAr 21; DoAr 21 may simply have moved into the vicinity of this gas because 
of its motion through the cloud, as we discuss further below. The FUV flux from DoAr 21 
has created a small-scale PDR in t his nearby gas, ex citing the observed H 2 , PAH, and VSG 
emission. As in models of PDRs (jAUers et al.ll2005l ). the 1-0 S(l) emis sion is at the edge 
of the PDR and extends farther into the ionized region, explaining why IBary et al.l (120021 ) 
det ect some emission i n an E-W slit that does not intersect the bright ridge to the north, 
but Bitner et al. ( 2008 ) do not detect emission in a transition that traces cooler gasl"1 



11 While the bright spot l'.'l to the N lies outside the slit used bv lBarv et al.l (|2002l ). if we assume a Gaussian 
profile for the seeing with 1'.'4 FWHM, then roughly 12% of the flux of a point source l'.'l away from the 
star would fall within the slit, so the northern ridge could still make some contribution to the observed flux. 
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One major testable prediction of this scenario is the spatial distribution of H2 emis- 
sion. This could be mapped with narrow-band filters, or (ideally) with an integral-field 
spectrograph that would also measure the line widths and line ratios to probe the physical 
conditions in the emitting gas. For example, observations of both the S(l) 1-0 and the S(l) 
2-1 lines would probe whether the gas is predominantly thermally excited or pumped by 
FUV photons. 

Spatially-unresolved observations at high spectral resolution should detect more flux 
in a N— S slit than in an E-W one; making observations close in time with different slit 
orientations would help determine whether time variability or spatial distribution is more 
important. Given that PDR models predict that 1-0 S(l) emission is spread over a range 
of emitting temperatures, the 1-0 S(l) line measured in a N -S slit may have a stronger 
flux and a somewhat narrower width than the line observed by iBary et al.l (120021 ) as it may 
sample both hotter and cooler gas. In contrast, if there is a Keplerian disk that produces 
the 1-0 S(l) emission, as suggested by IBary et al.l (120021 ). then closely-spaced observations 
with different slit orientations should produce the same flux and line width, and mapping of 
the line emission should show it closely associated with the stellar position of DoAr 21. 

The proper motion of DoAr 21 provides an interesting, though admittedly speculative, 
piece of circumstantial evidence in favor of the hypothes is that DoAr 2 1 is i lluminating; 
ambient cloud material rather than a circumstellar disk. lLoinard et al.l (120081 ) measure a 
proper motion of (/i a cos5,n s ) = (-26.47 ± 0.92,-28.23 ± 0.73) mas yr" 1 for DoAr 21. 
This proper motion in RA is significantly different from th e mean proper motion for p Oph 
members of ( — 10 ± 2, —27 ± 2) mas yr -1 (IMamajeki 120081 ). and yet there is ample evidence 
from the measured distance, luminosity, Li absorption, and X -ray properties that DoAr 21 
is a young star associated with the cloud. iMamajeki (120081 ) gives proper motions for 38 
candidate p Oph members. Not counting the three stars that Mamajek rejects as non- 
members, the only star with an RA proper motion as disc repant from the mean as DoAr 
21's motion is GSS 20 (ROX 7), which Iwilking et al.l (b()05h flag as "dwarf?", i.e. a possible 
non- member. 

If the mean p Oph proper motion is representative of the space motion of the gas near 
DoAr 21, then DoAr 21 has a tangential motion of 16.5 ± 2.2 mas yr -1 (9.5 ± 1.3 km s _1 ) 
westward relative to the gas. This is consistent with the fact that there is less material 
visible to the east of DoAr 21 (Fig. H]) than in other directions. At this velocity, DoAr 21 
would have been 1" to the east just 60 years ago, and could have photoevaporated or ionized 
much of the material in its vicinity, just as it appears to be doing now in the region within 
100 AU of the star. Such a high velocity relative to the ambient gas would be unusual, given 
the typical 1-2 km s" 1 velocity dispersions of young associations, but not unheard of. One 
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example is the high- velocity young star PV Cep (jGoodman fe Arceil2004l ). which is moving 
at ~ 20 km s _1 relative to its local cloud. 

To close, we note that the strong PAH emission at large distances from DoAr 21, and 
the implied strength of the FUV flux needed to excite it, lends empirical support to the hy- 
pothesis of I Alexander. Clarke, fc Pringld (120051 ) that there may be significant FUV flux from 
T Tauri stars even after accretion has ceased. Such a flux might enable photoevaporation to 
disperse outer disks on a short timescale even after accretion from the inner disk ceases. 



5. Conclusions 



We have presented a variety of new data on DoAr 21 in an attempt to understand its 
circumstellar environment and the influence of high-energy radiation on it. High- resolution 
mid-infrared images show PAH and very small grain emission in an irregular distribution 
over hundreds of AU from the star, with little emission from the inner 100 AU where a 
circumstellar disk might lie. A new high-resolution X-ray grating spectrum from Chandra 
reveals a very luminous, hard corona, dominated by plasma at temperatures up to 10 8 K; 
neither this spectrum nor new high-resolution optical spectra show any evidence for accretion. 
The inferred FUV radiation from the transition region of DoAr 21 is strong enough to excite 
the observed PAH, VSG, and H 2 emission, creating a small-scale PDR around DoAr 21. 
This strong radiation may also have played a role in removing any disk from around DoAr 
21, despite its young age. 
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